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We review the general relativistic theory of the motion, and of the timing, 
of binary systems containing compact objects (neutron stars or black holes). 
Then we indicate the various ways one can use binary pulsar data to test the 
strong-field and/or radiative aspects of General Relativity, and of general classes 
of alternative theories of relativistic gravity. 

1 Introduction 

The discovery of binary pulsars in 1974 [T] opened up a new testing ground 
for relativistic gravity. Before this discovery, the only available testing ground 
for relativistic gravity was the solar system. As Einstein's theory of General 
Relativity (GR) is one of the basic pillars of modern science, it deserves to 
be tested, with the highest possible accuracy, in all its aspects. In the solar 
system, the gravitational field is slowly varying and represents only a very small 
deformation of a flat spacetime. As a consequence, solar system tests can only 
probe the quasi-stationary (non radiative) weak-field limit of relativistic gravity. 
By contrast binary systems containing compact objects (neutron stars or black 
holes) involve spacetime domains (inside and near the compact objects) where 
the gravitational field is strong. Indeed, the surface relativistic gravitational 
field Hqq ~ 2 GM/c^R of a neutron star is of order 0.4, which is close to the one 
of a black hole (2 GM/c^R = 1) and much larger than the surface gravitational 
fields of solar system bodies: (2 GM/c'^R)sun - 10"^ (2 GM/c^ i?)Earth - 10"^. 
In addition, the high stability of "pulsar clocks" has made it possible to monitor 
the dynamics of its orbital motion down to a precision allowing one to measure 
the small (~ {v/c)^) orbital effects linked to the propagation of the gravitational 
field at the velocity of light between the pulsar and its companion. 

The recent discovery of the remarkable double binary pulsar PSR J0737— 
3039 [21 [3] (see also the contributions of M. Kramer and A. Possenti to these 
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proceedings) has renewed the interest in the use of binary pulsars as test-beds 
of gravity theories. The aim of these notes is to provide an introduction to the 
theoretical frameworks needed for interpreting binary pulsar data as tests of GR 
and alternative gravity theories. 

2 Motion of binary pulsars in general relativity 

The traditional (text book) approach to the problem of motion of N separate 
bodies in GR consists of solving, by successive approximations, Einstein's field 
equations (we use the signature — h ++) 

1 _ SttG 

together with their consequence 

T^"^ = . (2) 
To do so, one assumes some specific matter model, say a perfect fluid, 

T^"' = {e+p)u^'u'' +pg^"' . (3) 
One expands (say in powers of Newton's constant) 

5m. (2^^ ) = '/m- + h^ilu + + . . . , (4) 

together with the use of the simplifications brought by the 'Post-Newtonian' 
approximation {dohf^u = c~^dthf^i, ^ dihf^u] v/c <^ 1, p <^ e). Then one 
integrates the local material equation of motion ([2|) over the volume of each 
separate body, labelled say by a = 1, 2, . . . , iV. In so doing, one must define 
some 'center of mass' of body a, as well as some (approximately conserved) 
'mass' vria of body a, together with some corresponding 'spin vector' 5* and, 
possibly, higher multipole moments. 

An important feature of this traditional method is to use a unique coor- 
dinate chart to describe the full A^-body system. For instance, the center 
of mass, shape and spin of each body a are all described within this common 
coordinate system x'^ . This use of a single chart has several inconvenient as- 
pects, even in the case of weakly self-gravitating bodies (as in the solar system 
case). Indeed, it means for instance that a body which is, say, spherically sym- 
metric in its own 'rest frame' X°' will appear as deformed into some kind of 
ellipsoid in the common coordinate chart x^^ . Moreover, it is not clear how 
to construct 'good definitions' of the center of mass, spin vector, and higher 
multipole moments of body a, when described in the common coordinate chart 
x^ . In addition, as we are interested in the motion of strongly self-gravitating 
bodies, it is not a priori justified to use a simple expansion of the type ([4]) be- 
cause h^i^iu ~ '}2Gma/{c^ \x — Za\) will not be uniformly small in the common 
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coordinate system . It will be small if one stays far away from each object a, 
but, as recalled above, it will become of order unity on the surface of a compact 
body. 

These two shortcomings of the traditional 'one-chart' approach to the rela- 
tivistic problem of motion can be cured by using a ^multi- chart' approach.The 
multi-chart approach describes the motion of N (possibly, but not necessarily, 
compact) bodies by using A^-|- 1 separate coordinate systems: (i) one global coor- 
dinate chart x'^ (/i = 0, 1, 2, 3) used to describe the spacetime outside N 'tubes', 
each containing one body, and (ii) N local coordinate charts (a = 0, 1, 2, 3; 
a — 1,2, .. . , N) used to describe the spacetime in and around each body a. 
The multi-chart approach was first used to discuss the motion of black holes 
and other compact objects [H [5l [6l [3 [H [9l [TOl |TT]. Then it was also found to 
be very convenient for describing, with the high-accuracy required for dealing 
with modern technologies such as VLBI, systems of iV weakly self-gravitating 
bodies, such as the solar system [T^ [T5] . 

The essential idea of the multi-chart approach is to combine the information 
contained in several expansions. One uses both a global expansion of the type 
(HI and several local expansions of the type 

G^0{X2) = G[3(X2; m,) + m„ m,) + • • • , (5) 

where Gl^0{X;ma) denotes the (possibly strong-field) metric generated by an 
isolated body of mass nia (possibly with the additional effect of spin). 

The separate expansions ^ and (O are then 'matched' in some overlapping 
domain of common validity of the type Gma/c^ < i?a <C \x — Za\ <C d ~ ja^a — a;b| 
(with b ^ a), where one can relate the different coordinate systems by expansions 
of the form 

x^ = z^{Ta) + ef (T,) + i /^(T,) X:Xi + --- (6) 

The multi-chart approach becomes simplified if one considers compact bodies 
(of radius Ra comparable to 2 Gnia/c'^). In this case, it was shown i9 , by 
considering how the 'internal expansion' ([5]) propagates into the 'external' one 
([4| via the matching ([6]), that, in General Relativity, the internal structure of 
each compact body was effaced to a very high degree, when seen in the external 
expansion ((4]). For instance, for non spinning bodies, the internal structure of 
each body (notably the way it responds to an external tidal excitation) shows 
up in the external problem of motion only at the fifth post-Newtonian (5PN) 
approximation, i.e. in terms of order (w/c)^*^ in the equations of motion. 

This 'effacement of internal structure ' indicates that it should be possible 
to simplify the rigorous multi-chart approach by skeletonizing each compact 
body by means of some delta-function source. Mathematically, the use of dis- 
tributional sources is delicate in a nonlinear theory such as GR. However, it 
was found that one can reproduce the results of the more rigorous matched- 
multi-chart approach by treating the divergent integrals generated by the use 
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of delta- function sources by means of (complex) analytic continuation [5] . The 
most efficient method (especially to high PN orders) has been found to use 
analytic continuation in the dimension of space d [14] . 

Finally, the most efficient way to derive the general relativistic equations of 
motion of N compact bodies consists of solving the equations derived from the 
action (where g = — det((7^i,)) 

s = J ^ Va Y^:^ Ria) - X! ™° ^ / \l -at^-i^iza) dz'i dz^ , (7) 

formally using the standard weak- field expansion ^ , but considering the space 
dimension d as an arbitrary complex number which is sent to its physical value 
d = 3 only at the end of the calculation. 

Using this methoc(3 one has derived the equations of motion of two compact 
bodies at the 2.5PN (w^/c^) approximation level needed for describing binary 
pulsars [UlIlSlE]: 

= Al„{Za- Zb)+C-^ Al^{Za- Zb,Va,Vb) 
+ C^"^ Al^4^{Za- Zb,Va,Vb,Sa,Sb) 

+ C-^Al,{Za~Zb,Va~Vb)+0{c~''). (8) 

Here AI^q — —Gmb{zl — zl)/\za — 2f,P denotes the Newtonian acceleration, A\2 
its IPN modification, A^^^ its 2PN modification (together with the spin-orbit 
effects), and the 2.5PN contribution of order v^/c^. [See the references 
above; or the review [17j . for more references and the explicit expressions of A2, 
A4 and A^.] It was verified that the term has the effect of decreasing the 
mechanical energy of the system by an amount equal (on average) to the energy 
lost in the form of gravitational wave flux at infinity. Note, however, that here 
was derived, in the near zone of the system, as a direct consequence of the 
general relativistic propagation of gravity, at the velocity c, between the two 
bodies. This highlights the fact that binary pulsar tests of the existence of 
are direct tests of the reality of gravitational radiation. 

Recently, the equations of motion ([8]) have been computed to even higher 
accuracy: 3PN ~ v^/c^ [HI [ISl EQI EH [22] and 3.5PN ~ uVc^ [23l Ell [25] 
(see also the review [2S])- These refinements are, however, not (yet) needed for 
interpreting binary pulsar data. 



3 Timing of binary pulsars in general relativity 

In order to extract observational effects from the equations of motion ([5]) one 
needs to go through two steps: (i) to solve the equations of motion ([8]) so as to 

^Or, more precisely, an essentially equivalent analytic continuation using the so-called 
'Riesz kernels'. 
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get the coordinate positions Z\ and Z2 as explicit functions of the coordinate 
time i, and (ii) to relate the coordinate motion Za(fy to the pulsar observables, 
i.e. mainly to the times of arrival of electromagnetic pulses on Earth. 

The first step has been accomplished, in a form particularly useful for dis- 
cussing pulsar timing, in Ref. [57]. There (see also [5S]) it was shown that, 
when considering the full (periodic and secular) effects of the A-i ^ jt? terms 
in Eq. ([5]), together with the secular effects of the A4 ~ jc^ and A5 ~ w^/c^ 
terms, the relativistic two-body motion could be written in a very simple 'quasi- 
Keplerian' form (in polar coordinates), namely: 

ndt ^ G ~ u — et sin u , (9) 



6* - 6*0 = (1 + fc) 2arctan 



u 

tan — 
2 



(10) 



R = ^-afc afl(l - e_RC0Su) , (11) 

Ta = \Za ~ Zcm\^ ar{\- erCO^u) , (12) 

Th = - ZcAf I = ar'(l - cosu) . (13) 

Here n = 2^/?}, denotes the orbital frequency, k = A0/27r = (tj)/n = 
{6j) Pb/27r the fractional periastron advance per orbit, u an auxiliary angle ('rel- 
ativistic eccentric anomaly'), et,ee,eR,er and e^' various 'relativistic eccentric- 
ities' and afi,ar and a^' some 'relativistic semi-major axes'. See [2 7) for the 
relations between these quantities, as well as their link to the relativistic energy 
and angular momentum E, J. A direct study [28j of the dynamical effect of the 
contribution A5 ~ u^/c^ in the equations of motion ([8]) has checked that it led 
to a secular increase of the orbital frequency n{t) ~ n{0) + n{t — to) , and thereby 
to a quadratic term in the 'relativistic mean anomaly' £ — J ndt + a appearing 
on the left-hand side (L.H.S.) of Eq. ©: 

i^ao + no{t~to) + ^n{t-to)^ . (14) 

As for the contribution A4 ~ w'*/c^ it induces several secular effects in the 
orbital motion: various 2PN contributions to the dimensionless periastron pa- 
rameter k ((54 fc ^ v'^/c^+ spin-orbit effects), and secular variations in the incli- 
nation of the orbital plane (due to spin-orbit effects). 

The second step in relating ^ to pulsar observations has been accomplished 
through the derivation of a 'relativistic timing formula ' [291 130j . The 'timing 
formula' of a binary pulsar is a multi-parameter mathematical function relating 
the observed time of arrival (at the radio-telescope) of the center of the N^^ 
pulse to the integer N. It involves many different physical effects: (i) dispersion 
effects, (ii) travel time across the solar system, (iii) gravitational delay due to 
the Sun and the planets, (iv) time dilation effects between the time measured 
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on the Earth and the solar-system-barycenter time, (v) variations in the travel 
time between the binary pulsar and the solar-system barycenter (due to relative 
accelerations, parallax and proper motion), (vi) time delays happening within 
the binary system. We shall focus here on the time delays which take place 
within the binary system (see the lectures of M. Kramer for a discussion of the 
other effects). 

For a proper derivation of the time delays occurring within the binary sys- 
tem we need to use the multi-chart approach mentionned above. In the 'rest 
frame' {X^ = cTa, X^) attached to the pulsar a, the pulsar phenomenon can be 
modelled by the secularly changing rotation of a beam of radio waves: 

= y" ^aiTa) dTa ~ ^^Ta + ^^tlaT^ + ^^aT^ + ■ ■ ■ , (15) 

where $a is the longitude around the spin axis. [Depending on the precise defi- 
nition of the rest-frame attached to the pulsar, the spin axis can either be fixed, 
or be slowly evolving, see e.g. [13 .] One must then relate the initial direction 
(0a, $a), and proper time T^, of emission of the pulsar beam to the coordinate 
direction and coordinate time of the null geodesic representing the electromag- 
netic beam in the 'global' coordinates x^^ used to describe the dynamics of the 
binary system [NB: the explicit orbital motion ((9|l- (fT3|) refers to such global 
coordinates x*^ = ci, x*]. This is done by using the link ([6]) in which denotes 
the global coordinates of the 'center of mass' of the pulsar, Ta the local (proper) 
time of the pulsar frame, and where, for instance 



°-^-i(l + i^ + 3^ + . ..)+... (16) 

c V 2 c^rab ' 



Using the link ^ (with expressions such as for the coefficients ef , . . .) 
one finds, among other results, that a radio beam emitted in the proper direction 
iV' in the local frame appears to propagate, in the global frame, in the coordinate 
direction where 

n' = N' + -- —— + 0[%] . (17) 

c c V c / 

This is the well known 'aberration effect', which will then contribute to the 
timing formula. 

One must also write the link between the pulsar 'proper time' Tq and the 
coordinate time t = x°/c — z^/c used in the orbital motion ([9l)- p3|) . This reads 

-^dTl^~g^,{al)dzidz'^^ (18) 

where the 'tilde' denotes the operation consisting (in the matching approach) in 
discarding in the 'self contributions' ^ {Gnia/ Ra)"' , while keeping the effect 
of the companion (^ Gnib/rab, etc. . .). One checks that this is equivalent (in 
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the dimensional-continuation approach) in taking x'^ = zf^ for sufficiently small 
values of the real part of the dimension d. To lowest order this yields the link 

which combines the special relativistic and general relativistic time dilation 
effects. Hence, following j30] we can refer to them as the 'Einstein time delay'. 

Then, one must compute the (global) time taken by a light beam emitted 
by the pulsar, at the proper time Ta (linked to iomission by (fTO)) '). in the initial 
global direction (see Eq. ([IT])), to reach the barycenter of the solar system. 
This is done by writing that this light beam follows a null geodesic: in particular 

= ds^ = g^A^^) dxf" dx" dt^ + (l + ^) (20) 

where U — Gma/\x — Za\+ Grni,/\x — zii\ is the Newtonian potential within the 
binary system. This yields (with = <cmission, ia = tanivai) 

ta-te= r dtc,^ \dx\ + 4 f + r^^) M^l • (21) 

Jt, ^ Jt, \\X-Za\ \X-Zb\j 

The first term on the last RHS of Eq. (pij) is the usual 'light crossing time' 
\ |-2barycenter(ia) " Za{te)\ between the pulsar and the solar barycenter. It con- 
tains the 'Roemer time delay' due to the fact that Za{te) moves on an orbit. 
The second term on the last RHS of Eq. (pij) is the 'Shapiro time delay' due to 
the propagation of the beam in a curved spacetime (only the Gmi, piece linked 
to the companion is variable). 

When inserting the 'quasi-Keplerian' form ((9|)- (fT3|) of the relativistic motion 
in the 'Roemer' term in (EJ), together with all other relativistic effects, one finds 
that the final expression for the relativistic timing formula can be significantly 
simplified by doing two mathematical transformations. One can redefine the 
'time eccentricity' appearing in the 'Kepler equation' ([5]), and one can define 
a new 'eccentric anomaly' angle: u — *■ [we henceforth drop the superscript 
'new' on u]. After these changes, the binary-system part of the general relativis- 
tic timing formula [30] takes the form (we suppress the index a on the pulsar 
proper time Ta) 

tbaryccnter " ^0 = ^"M?^ + Afl(T) + Ae{T) + As{T) + Aa{T)] (22) 

with 

— xsin(jj[cosu — e(l -I- (5j.)] + a;[l — e^(l -|- (Je)^]^^^ coswsinu , (23) 
Ae = 7sinM, (24) 
A5 — — 2rln{l — ecosu — s[sinci;(cosM — e) -I- (1 — e^)"'^/^ coswsinM]}(25) 
Aa = A{s'm[LLj + Ae{u)] + e sin uj} + B {cos[uj + Ae{u)] + e cos Lu} , (26) 
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where x = xo + x(T — Tq) represents the projected hght-crossing time {x = 
Opuisar sini/c), 6 = cq + e{T — Tq) a certain (relativistically-defined) 'timing 
eccentricity', Ae{u) the function 



Ae{u) = 2 arctan 



tan — 



1 -e / 2 



(27) 



Lo — loq + k Ae{u) the 'argument of the periastron', and where the (relativisti- 
cally-defined) 'eccentric anomaly' u is the function of the 'pulsar proper time' 
T obtained by solving the Kepler equation 



e sin u ^ 2tt 



Pb 2 ' I Pfc 



(28) 



It is understood here that the pulsar proper time T corresponding to the N^^ 
pulse is related to the integer N by an equation of the form 

N ^co + VpT +]^UpT^ ^^DpT"^ . (29) 

From these formulas, one sees that 5e (and 5r) measure some relativistic distor- 
tion of the pulsar orbit, 7 the amplitude of the 'Einstein time delayl As, and 
r and s the range and shape of the 'Shapiro time delayH A5. Note also that 
the dimensionless PPK parameter k measures the non-uniform advance of the 
periastron. It is related to the often quoted secular rate of periastron advance 
u! = (duj/dt) by the relation k = loP^ /2tt. It has been explicitly checked that 
binary-pulsar observational data do indeed require to model the relativistic pe- 
riastron advance by means of the non-uniform (and non-trivial) function of u 
multiplying k on the R.H.S. of Eq. ^ ^310. Finally, we see from Eq. ([28]) that 
Pb represents the (periastron to periastron) orbital period at the fiducial epoch 
To, while the dimensionless parameter Pf, represents the time derivative of Ph 
(at To). 

Schematically, the structure of the DD timing formula ([22|) is 

Wontcr -h = F [Tat; {/^'}; {p^^}; {g^^}] , (30) 

where ibaryccntor dcnotcs the solar-system barycentric (infinite frequency) ar- 
rival time of a pulse, T the pulsar emission proper time (corrected for aberra- 
tion), {p^} = {Pb, To, Co, Wo, xq} is the set of Keplerian parameters, {p^^ = 
k, 7, Pb, r, s, Sg, e, x} the set of separately measurable post-Keplerian parameters. 



^Thc post-Keplerian timing parameter 7, first introduced in 1291 . has tiie dimension of 
time, and should not be confused with the dimensionless post-Newtonian Eddington parameter 
^PPN pj-Qbed by solar-system experiments (see below). 

^The dimensionless parameter s is numerically equal to the sine of the inclination angle i 
of the orbital plane, but its real definition within the PPK formalism is the timing parameter 
which determines the 'shape' of the logarithmic time delay As(T). 

''Alas this function is theory-independent, so that the non-uniform aspect of the periastron 
advance cannot be used to yield discriminating tests of relativistic gravity theories. 
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and {q^^} = {5r, A, B, D} the set of not separately measurable post-Keplerian 
parameters [31]. [The parameter D is a 'Doppler factor' which enters as an 
overall multiplicative factor on the right-hand side of Eq. ([H]).] 

A further simplification of the DD timing formula was found possible. In- 
deed, the fact that the parameters {q^^} = {5r,A,B,D} are not separately 
measurable means that they can be absorbed in changes of the other param- 
eters. The explicit formulas for doing that were given in [30] and [31]: they 
consist in redefining e,x,Pb,Se and 6r. At the end of the day, it suffices to 
consider a simplified timing formula where {Sr, A, B, D} have been set to some 
given fiducial values, e.g. {0, 0, 0, 1}, and where one only fits for the remaining 
parameters {p^} and {p^^}. 

Finally, let us mention that it is possible to extend the general parametrized 
timing formula by writing a similar parametrized formula describing the ef- 
fect of the pulsar orbital motion on the directional spectral luminosity [d(energy) 
/c?(time) d(frequency) rf(solid angle)] received by an observer. As discussed in 
detail in [21] this introduces a new set of 'pulse-structure post-Keplerian pa- 
rameters'. 

4 Phenomenological approach to testing rela- 
tivistic gravity with binary pulsar data 

As said in the Introduction, binary pulsars contain strong gravity domains and 
should therefore allow one to test the strong-field aspects of relativistic gravity. 
The question we face is then the following: How can one use binary pulsar data 
to test strong-field (and radiative) gravity? 

Two different types of answers can be given to this question: a phenomeno- 
logical (or theory-independent) one, or various types of theory- dependent ap- 
proaches. In this Section we shall consider the phenomenological approach. 

The phenomenological approach to binary-pulsar tests of relativistic gravity 
is called the parametrized post-Keplerian formalism [321 I31| . This approach 
is based on the fact that the mathematical form of the multi-parameter DD 
timing formula (|30p was found to be applicable not only in General Relativity, 
but also in a wide class of alternative theories of gravity. Indeed, any theory 
in which gravity is mediated not only by a metric field 17^,^ but by a general 
combination of a metric field and of one or several scalar fields ip^""^ will induce 
relativistic timing effects in binary pulsars which can still be parametrized by the 
formulas ((22|) - (|29p . Such general 'tensor- multi-scalar' theories of gravity contain 
arbitrary functions of the scalar fields. They have been studied in full generality 
in [33] . It was shown that, under certain conditions, such tensor-scalar gravity 
theories could lead, because of strong-field effects, to very different predictions 
from those of General Relativity in binary pulsar timing observations [Ml [551 [55] . 
However, the point which is important for this Section, is that even when such 
strong-field effects develop one can still use the universal DD timing formula 
([30]) to fit the observed pulsar times of arrival. 
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The basic idea of the phenomenological, parametrized post-Keplerian (PPK) 
approach is then the following: By least-square fitting the observed sequence 
of pulsar arrival times to the parametrized formula (|30p (in which T/v is 
defined by Eq. (P^)) which introduces the further parameters Vp, Vp, Up) one can 
phenomenologically extract from raw observational data the (best fit) values of 
all the parameters entering Eqs. ([M)) and (|30p . In particular, one so determines 
both the set of Keplerian parameters {p-^^ } = {Pf,, Tq, cq, wq, xq}, and the set of 
post-Keplerian (PK) parameters {p^^} = {k,j,Pb,r,s,Se,e,±}. In extracting 
these values, we did not have to assume any theory of gravity. However, each 
specific theory of gravity will make specific predictions relating the PK param- 
eters to the Keplerian ones, and to the two (a priori unknown) masses TOq and 
rrib of the pulsar and its companion. [For certain PK parameters one must also 
consider other variables related to the spin vectors of a and b.] In other words, 
the measurement (in addition of the Keplerian parameters) of each PK param- 
eter defines, for each given theory, a curve in the {ma,mh) mass plane. For any 
given theory, the measurement of two PK parameters determines two curves 
and thereby generically determines the values of the two masses nia and mf, (as 
the point of intersection of these two curves). Therefore, as soon as one mea- 
sures three PK parameters one obtains a test of the considered gravity theory. 
The test is passed only if the three curves meet at one point. More generally, 
the measurement of n PK timing parameters yields n — 2 independent tests of 
relativistic gravity. Any one of these tests, i.e. any simultaneous measurement 
of three PK parameters can either confirm or put in doubt any given theory of 
gravity. 

As General Relativity is our current most successful theory of gravity, it is 
clearly the prime target for these tests. Wc have seen above that the timing data 
of each binary pulsar provides a maximum of 8 PK parameters: k, 7, Pb, r, s, Se, e 
and X. Here, we were talking about a normal 'single line' binary pulsar where, 
among the two compact objects a and b only one of the two, say a is observed 
as a pulsar. In this case, one binary system can provide up to 8 — 2 = 6 tests of 
GR. In practice, however, it has not yet been possible to measure the parameter 
dg (which measures a small relativistic deformation of the elliptical orbit), nor 
the secular parameters e and x. The original Hulse- Taylor system PSR 1913-1-16 
has allowed one to measure 3 PK parameters: k = (cO)Pb/27r, 7 and Pb- The 
two parameters k and 7 involve (non radiative) strong-field effects, while, as 
explained above, the orbital period derivative Pb is a direct consequence of the 
term A5 ~ w^/c^ in the binary-system equations of motion ([5]). The term ^5 is 
itself directly linked to the retarded propagation, at the velocity of light, of the 
gravitational interaction between the two strongly self-gravitating bodies a and 
b. Therefore, any test involving Pb will be a mixed radiative strong-field test. 

Let us explain on this example what information one needs to implement a 
phenomenological test such as the (A:— 7— f';,)i9i3+i6 one. First, we need to know 
the predictions made by the considered target theory for the PK parameters fc, 7 
and Pb as functions of the two masses and rub- These predictions have been 
worked out, for General Relativity, in Refs. [29l[28l[30]. Introducing the notation 
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(where n = 2-k / Pb) 

M = ma + nib (31) 

Xa = nia/M ; Xb = rub/M ; X, + Xb=l (32) 

. ,33) 



they read 



fc^^(m„m,,) = -2/3^, (34) 



e 



7^^^(m„mb) - - Xb(l + X6)/3^, (35) 



n 



pGR/ \ _ 19^^ 1 + 24 + 96 6^ \- \- «5 /qf?^ 
^6 [ma,mb) (1 - e^)"^/^ XaXbPo- (oo) 

However, if we use the three predictions (p4|) - ((36)) . together with the best 
current observed values of the PK parameters /jobs^ ^obs^ pobd ^^y^ shall find 

that the three curves k'^^{ma, mb) = j'^^{nia,mb) = 7°'^^ P^^{ma, mb) = 
P^^^ in the {nia^mb) mass plane fail to meet at about the 13 ct level! Should 
this put in doubt General Relativity? No, because Ref. [38] has shown that the 
time variation (notably due to galactic acceleration effects) of the Doppler fac- 
tor D entering Eq. entailed an extra contribution to the 'observed' period 
derivative P^^^ . We need to subtract this non-GR contribution before drawing 
the corresponding curve: P^^{ma,mb) = Pb^^ — pgaiactic^ Then one finds that 
the three curves do meet within one cr. This yields a deep confirmation of Gen- 
eral Relativity, and a direct observational proof of the reality of gravitational 
radiation. 

We said several times that this test is also a probe of the strong-field aspects 
of GR. How can one see this? A look at the GR predictions p4l) - (p6| does not 
exhibit explicit strong-field effects. Indeed, the derivation of Eqs. used 
in a crucial way the 'effacement of internal structure' that occurs in the general 
relativistic dynamics of compact objects. This non trivial property is rather 
specific of GR and means that, in this theory, all the strong- field effects can be 
absorbed in the definition of the masses ma and mb. One can, however, verify 
that strong-field effects do enter the observable PK parameters fc, 7, Pb etc. . . by 
considering how the theoretical predictions get modified in alternative 

theories of gravity. The presence of such strong-field effects in PK parameters 
was first pointed out in Ref. [7] (see also [S^) for the Jordan-Fierz-Brans-Dicke 
theory of gravity, and in Ref. [5] for Rosen's bi-metric theory of gravity. A 
detailed study of such strong-field deviations was then performed in [33, 341135] 
for general tensor-(multi-)scalar theories of gravity. In the following Section 
we shall exhibit how such strong-field effects enter the various post-Keplerian 
parameters. 

Continuing our historical review of phcnomenological pulsar tests, let us 
come to the binary system which was the first one to provide several 'pure 



11 



strong-field tests' of relativistic gravity, without mixing of radiative effects: 
PSR 1534+12. In this system, it was possible to measure the four (non ra- 
diative) PK parameters k,j,r and s. [We see from Eq. (|25p that r and s mea- 
sure, respectively, the range and the shape of the 'Shapiro time delay' As-] The 
measurement of the 4 PK parameters k,'y,r,s define 4 curves in the {ma,rnb) 
mass plane, and thereby yield 2 strong-field tests of GR. It was found in [JU] 
that GR passes these two tests. For instance, the ratio between the measured 
value 5°'^'' of the phenomenological parameter s and the value s^^[/e°'^'', 7°'^''] 
predicted by GR on the basis of the measurements of the two PK parameters k 
and 7 (which determine, via Eqs. ([M)) , the GR-predicted value of ma and 
mb) was found to be s°''Vs^^ [/^°''^ 7°*^"] = 1.004 ± 0.007 00]. The most recent 
data [H] yield s°''Vs°^ [/c°''^ 7°'""] = 1.000 ± 0.007. We see that we have here a 
confirmation of the strong-field regime of GR at the 1% level. 

Another way to get phenomenological tests of the strong field aspects of 
gravity concerns the possibility of a violation of the strong equivalence principle. 
This is parametrized by phenomenologically assuming that the ratio between the 
gravitational and the inertial mass of the pulsar differs from unity (which is its 
value in GR): ('Tigi.av/winGrt)a = Similarly to what happens in the Earth- 

Moon-Sun system [32], the three-body system made of a binary pulsar and of the 
Galaxy exhibits a 'polarization' of the orbit which is proportional to A = — 
Af), and which can be constrained by considering certain quasi-circular neutron- 
star- white-dwarf binary systems [33] • See [33] for recently published improved 
hmitfS on the phenomenological equivalence-principle violation parameter A. 

The Parkes multibeam survey has recently discovered several new interesting 
'relativistic' binary pulsars, thereby giving a huge increase in the number of 
phenomenological tests of relativistic gravity. Among those new binary pulsar 
systems, two stand out as superb testing grounds for relativistic gravity: (i) 
PSR J1141-6545 [461 Hz], and (ii) the remarkable double binary pulsar PSR 
J0737-3039A and B [1 [31 [31133] (see also the lectures by M. Kramer and 
A. Possenti). 

The PSR J1141— 6545 timing data have led to the measurement of 3 PK 
parameters: k, 7, and Pb [47]. As in PSR 1913+16 this yields one mixed 
radiative-strong-field tes10. 



^As already mentioned the dimensionless parameter s is numerically equal (in all theories) 
to the sine of the inclination angle i of the orbital plane, but it is better thought, in the PPK 
formalism, as a phenomenological timing parameter determining the 'shape' of the logarithmic 
time delay As(T). 

®Note, however, that these limits, as well as those previously obtained in |45) . assume 
that the (a priori pulsar-mass dependent) parameter A ~ Aa is the same for all the analyzed 
pulsars. 

'^In addition, scintillation data have led to an estimate of the sine of the orbital inclination, 
sini |50j . As said above, sinj numerically coincides with the PK parameter s measuring the 
'shape' of the Shapiro time delay. Therefore, one could use the scintillation measurements as 
an indirect determination of s, thereby obtaining two independent tests from PSR J1141— 6545 
data. A caveat, however, is that the extraction of sini from scintillation measurements rests 
on several simplifying assumptions whose validity is unclear. In fact, in the case of PSR 
J0737— 3039 the direct timing measurement of s disagrees with its estimate via scintillation 
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The timing data of the millisecond binary pulsar PSR J0737— 3039A have 
led to the direct measurement of 5 PK parameters: fc, 7, r, s and Pb [H |48l 
135] ■ In addition, the 'double line' nature of this binary system (i.e. the fact 
that one observes both components, A and B, as radio pulsars) allows one to 
perform new phenomenological tests by using Keplerian parameters. Indeed, the 
simultaneous measurement of the Keplerian parameters Xa and Xh representing 
the projected light crossing times of both pulsars {A and B) gives access to the 
combined Keplerian parameter 

^obs 

= ^ . (37) 

On the other hand, the general derivation of [3D] (applicable to any Lorentz- 
invariant theory of gravity, and notably to any tensor-scalar theory) shows that 
the theoretical prediction for the the ratio i?, considered as a function of the 
masses vria and TOf,, is 



i?'h™'T = + . (38) 



The absence of any explicit strong-field-gravity effects in the theoretical predic- 
tion ([55)1 (to be contrasted, for instance, with the predictions for PK parameters 
in tensor-scalar gravity discussed in the next Section) is mainly due to the con- 
vention used in |30| and ^31j for defining the masses and m;,. These are 
always defined so that the Lagrangian for two non interacting compact objects 
reads — nia c^(l — t^^/c^)^/^. In other words, nia <? represents the total 

a 

energy of body a. This means that one has implicitly lumped in the definition 
of rria many strong-self-gravity effects. [For instance, in tensor-scalar gravity 
rria includes not only the usual Einsteinian gravitational binding energy due 
to the self-gravitational field g^^{x), but also the extra binding energy linked 
to the scalar field f{x).] Anyway, what is important is that, when performing 
a phenomenological test from the measurement of a triplet of parameters, e.g. 
{fc,7,i?}, at least one parameter among them be a priori sensitive to strong- 
field effects. This is enough for guaranteeing that the crossing of the three 
curves k^^''"'^ {ma^mb) = A;°^% 7"^™'^y(TOa, Wfc) = 7°''", i?*''™'^ (m^, mfc) = R"^" 
is really a probe of strong-field gravity. 

In conclusion, the two recently discovered binary pulsars PSR J1141— 6545 
and PSR J0737— 3039 have more than doubled the number of phenomenological 
tests of (radiative and) strong-field gravity. Before their discovery, the 'canoni- 
cal' relativistic binary pulsars PSR 1913+16 and PSR 1534-1-12 had given us four 



data 1491 . It is therefore safer not to use scintillation estimates of sin i on the same footing as 
direct timing measurements of the PK parameter s. On the other hand, a safe way of obtaining 
an s-related gravity test consists in using the necessary mathematical fact that s = sini < 1. 
In GR the definition Xa = Qa sini/c leads to sini = nXa/iPo Xj). Therefore we can write the 
inequality nXa/{f3o{M) Xb) < 1 as a phenomenological test of GR. 
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such tests: one (fc — 7 — Pfc) test from PSR 1913+16 and three (fc — 7 — r — s — pJl) 
tests from PSR 1534+12. The two new binary systems have given us more 
phenomenological tests: one {k — j — Pb) (or two, k — j — Pb — s) tests from PSR 
J1141-6545 and four (fc - 7 - r - s - A - -R) tests from PSR J0737-303d3. As 
illustrated in Figure [U these nine phenomenological tests of strong-field (and 
radiative) gravity are all in beautiful agreement with General Relativity. 

In addition, let us recall that several quasi-circular wide binaries, made of 
a neutron star and a white dwarf, have led to high-precision phenomenological 
confirmations [H] (in strong-field conditions) of one of the deep predictions of 
General Relativity: the 'strong' equivalence principle, i.e. the fact that var- 
ious bodies fall with the same acceleration in an external gravitational field, 
independently of the strength of their self-gravity. 

Finally, let us mention that Ref. [31] has extended the philosophy of the 



^The timing measurement of P^^" in PSR 1534+12 is even more strongly affected by 
kinematic corrections (D terms) than in the PSR 1913+16 case. In absence of a precise, 
independent measurement of the distance to PSR 1534+12, the fc — 7 — Pj, test yields, at best, 
a ~ 15% test of GR. 

^Or even six, if we use the scintillation determination of s in PSR J1141— 6545. 
^"The companion pulsar 0737— 3039B being non recycled, and being visible only during a 
small part of its orbit, cannot be timed with sufficient accuracy to allow one to measure any 
of its post-Keplerian parameters. 
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phenomenological (parametrized post-Keplerian) analysis of timing data, to a 
similar phenomenological analysis of pulse- structure data. Ref. |31| showed that, 
in principle, one could extract up to 11 'post-Keplerian pulse-structure param- 
eters'. Together with the 8 post-Keplerian timing parameters of a (single-line) 
binary pulsar, this makes a total of 19 phenomenological PK parameters. As 
these parameters depend not only on the two masses rria, rUb but also on the two 
angles A, rj determining the direction of the spin axis of the pulsar, the maximum 
number of tests one might hope to extract from one (single-line) binary pulsar 
is 19 — 4 15. However, the present accuracy with which one can model and 
measure the pulse structure of the known pulsars has not yet allowed one to 
measure any of these new pulse-structure parameters in a theory-independent 
and model-independent way. 

Nonetheless, it has been possible to confirm the reality (and order of mag- 
nitude) of the spin-orbit coupling in GR which was pointed out [521 153j to be 
observable via a secular change of the intensity profile of a pulsar signal. Confir- 
mations of general relativistic spin-orbit effects in the evolution of pulsar profiles 
were obained in several pulsars: PSR 1913-hl6 [51 [Hj, PSR B1534-I-12 [5^1 and 
PSR J1141— 6545 [57j. In this respect, let us mention that the spin-orbit interac- 
tion affects also several PK parameters, either by inducing a secular evolution in 
some of them (see |31]) or by contributing to their value. For instance, the spin- 
orbit interaction contributes to the observed value of the periastron advance 
parameter k an amount which is significant for the pulsars (such as 1913-1-16 
and 0737—3039) where k is measured with high-accuracy. It was then pointed 
out [58j that this gives, in principle, and indirect way of measuring the moment 
of inertia of neutron stars (a useful quantity for probing the equation of state 
of nuclear matter [59l I60j). However, this can be done only if one measures, 
besides k, two other PK parameters with 10~^ accuracy. A rather tall order 
which will be a challenge to meet. 

The phenomenological approach to pulsar tests has the advantage that it can 
confirm or invalidate a specific theory of gravity without making assumptions 
about other theories. Moreover, as General Relativity has no free parameters, 
any test of its predictions is a potentially lethal test. From this point of view, it is 
remarkable that GR has passed with flying colours all the pulsar tests if has been 
submitted to. [See, notably. Fig. [T]] As argued above, these tests have probed 
strong-field aspects of gravity which had not been probed by solar-system (or 
cosmological) tests. On the other hand, a disadvantage of the phenomenological 
tests is that they do not tell us in any precise way which strong- field structures, 
have been actually tested. For instance, let us imagine that one day one specific 
PPK test fails to be satisfied by GR, while the others are OK. This leaves us in a 
quandary: If we trust the problematic test, we must conclude that GR is wrong. 
However, the other tests say that GR is OK. This example shows that we would 
like to have some idea of what physical effects, linked to strong-field gravity, 
enter in each test, or even better in each PK parameter. The 'effacement of 
internal structure' which takes place in GR does not allow one to discuss this 
issue. This gives us a motivation for going beyond the phenomenological PPK 
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approach by considering theory-dependent formalisms in which one embeds GR 
within a space of alternative gravity theories. 

5 Theory-space approach to testing relativistic 
gravity with binary pulsar data 

A complementary approach to testing gravity with binary pulsar data consists 
in embedding General Relativity within a multi-parameter space of alternative 
theories of gravity. In other words, we want to contrast the predictions of GR 
with the predictions of continuous families of alternative theories. In so doing we 
hope to learn more about which structures of GR are actually being probed in 
binary pulsar tests. This is a bit similar to the well-known psycho-physiological 
fact that the best way to appreciate a nuance of colour is to surround a given 
patch of colour by other patches with slightly different colours. This makes it 
much easier to detect subtle differences in colour. In the same way, we hope to 
learn about the probing power of pulsar tests by seeing how the phenomeno- 
logical tests summarized in Fig. [1] fail (or continue) to be satisfied when one 
continuously deform, away from GR, the gravity theory which is being tested. 

Let us first recall the various ways in which this theory-space approach has 
been used in the context of the solar-system tests of relativistic gravity. 

5.1 Theory-space approaches to solar-system tests of rel- 
ativistic gravity 

In the quasi-stationary weak-field context of the solar-system, this theory-space 
approach has been implemented in two different ways. First, the parametrized 
post-Newtonian (PPN) formalism [SSillMlllllMliSliniES] describes 
many 'directions' in which generic alternative theories of gravity might dif- 
fer in their weak-field predictions from GR. In its most general versions the 
PPN formalism contains 10 'post-Einstein' PPN parameters, 7 = j^^^ — l ["1 . 
/3 = Z?^^'^ — 1 , ^, ai , q;2 , aa , Ci , C2 , Cs I C4 • Each one of these dimensionless quanti- 
ties parametrizes a certain class of slow-motion, weak-field gravitational effects 
which deviate from corresponding GR predictions. For instance, 7 parametrizes 
modifications both of the effect of a massive body (say, the Sun) on the light 
passing near it, and of the terms in the two-body gravitational Lagrangian which 
are proportional to (Gma rrib/rab) ■ (I'a — Vh)'^/c^. 

A second way of implementing the theory-space philosophy consists in con- 
sidering some explicit, parameter-dependent family of alternative relativistic 
theories of gravity. For instance, the simplest tensor-scalar theory of gravity 



^^The PPN parameter 7^^^ is usually denoted simply as 7. To distinguish it from the 
Einstein-time-delay PPK timing parameter 7 used above we add the superscript PPN. In 
addition, as the value of 7^'^''^ in GR is 1, we prefer to work with the parameter 7 = 7^^* — 1 
which vanishes in GR, and therefore measures a 'deviation' from GR in a certain 'direction' 
in theory-space. Similarly with /3 = /3^^^ — 1. 
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put forward by Jordan [ST], Fierz [BS] and Brans and Dicke [55] has a unique 
free parameter, say = (2wbd + 3)~^. When aQ — » 0, this theory reduces 
to GR, so that (o^ 1/^bd) measures all the deviations from GR. When 
considering the weak-field limit of the Jordan-Fierz-Brans-Dicke (JFBD) the- 
ory, one finds that it can be described within the PPN formalism by choosing 
7 = -2ag(l-hag)"\ /3 = and ^ = a* = = 0. 

Having briefly recalled the two types of theory-space approaches used to 
discuss solar-system tests, let us now consider the case of binary-pulsar tests. 

5.2 Theory-space approaches to binary-pulsar tests of rel- 
ativistic gravity 

There exist generalizations of these two different theory-space approaches to the 
context of strong-field gravity and binary pulsar tests. First, the PPN formalism 
has been (partially) extended beyond the 'first post-Newtonian' (IPN) order 
deviations from GR (~ w^/c^ -t- Gm/c^ r) to describe 2PN order deviations from 

GR (^p" + §t) ^ Remarkably, there appear only two new parameters 

at the 2PN leve{ll: e and C. Also, by expanding in powers of the self-gravity 
parameters of body a and b the predictions for the PPK timing parameters 
in generic tensor-multi-scalar theories, one has shown that these predictions 
depended on several 'layers' of new dimensionless parameters |33| . Early among 
these parameters one finds, the IPN parameters (3, 7 and then the basic 2PN 
parameters e and (, but one also finds further parameters P3, /?",... 
which would not enter usual 2PN effects. The two approaches that we have just 
mentionned can be viewed as generalizations of the PPN formalism. 

There exist also useful generalizations to the strong-field context of the idea 
of considering some explicit parameter-dependent family of alternative theo- 
ries of relativistic gravity. Early studies [3 [HI [35] focussed either on the one- 
parameter JFBD tensor-scalar theory, or on some theories which are not con- 
tinuously connected to GR, such as Rosen's bimetric theory of gravity. Though 
the JFBD theory exhibits a marked difference from GR in that it predicts the 
existence of dipole radiation, it has the disadvantage that the weak field, solar- 
system constraints on its unique parameter are so strong that they drastically 
constrain (and essentially forbid) the presence of any non-radiative, strong-field 
deviations from GR. In view of this, it is useful to consider other 'mini-spaces' 
of alternative theories. 

A two-parameter mini-space of theories, that we shall denot^^ here as 
T2{P' , P"), was introduced in [33]. This two-parameter family of tensor-bi-scalar 
theories was constructed so as to have exactly the same first post-Newtonian 
limit as GR (i.e. 7 = /3 = • • • = 0), but to differ from GR in its predictions 

^^When restricting oneself to the general class of tensor-multi-scalar theories. At the IPN 
level, this restriction would imply that only the 'directions' 7 and /3 are allowed. 

^^We add here an index 2 to T as a reminder that this is a class of tensor- 6j-scalar theories, 
i.e. that they contain two independent scalar fields 1^1,1^2 besides a dynamical metric gpi/. 
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for the various observables that can be extracted from binary pulsar data. Let 
us give one example of this behaviour of the T2{f3',f3") class of theories. For 
a general theory of gravity we expect to have violations of the strong equiva- 
lence principle in the sense that the ratio between the gravitational mass of a 
self-gravitating body to its inertial mass will admit an expansion of the type 

^grav 2 

= l + Aa^l--mCa+V2cl + ... (39) 

where Ca = —2 ^J"^™" measures the 'gravitational compactness' (or fractional 
gravitational binding energy, Ca — — 2 E'f'^^/ma c^) of body a. The numerical 
coefficient rji of the contribution linear in Ca is a combination of the first post- 
Newtonian order PPN parameters, namely 771 = 4/3 — 7 [42 . The numerical 
coefficient 772 of the term quadratic in Ca is a combination of the IPN and 2PN 
parameters. When working in the context of the T2{f3',P") theories, the IPN 
parameters vanish exactly {P — — ^) and the coefficient of the quadratic term 
becomes simply proportional to the theory parameter fi' '■ ri2 = \BP\ where 
B « 1.026. This example shows explicitly how binary pulsar data (here the data 
constraining the equivalence principle violation parameter A — — At,, see 
above) can go beyond solar-system experiments in probing certain strong-self- 
gravity effects. Indeed, solar-system experiments are totally insensitive to 2PN 
parameters because of the smallness of Ca ~ Giria/c^ Ra and of the structure 
of 2PN effects [70]. By contrast, the 'compactness' of neutron stars is of order 
Ca ~ 0.21 rria/MQ ~ 0.3 [33j so that the pulsar limit | A| < 5.5 x 10"^ M} yields, 
within the T2{f3',f3") framework, a significant limit on the dimensionless (2PN 
order) parameter /?' : < 0.12. 

Ref. [35] introduced a new two-parameter mini-space of gravity theories, 
denoted here as Ti(ao7 Po), which, from the point of view of theoretical physics, 
has several advantages over the T2{(3',(3") mini-space mentionned above. First, 
it is technically simpler in that it contains only one scalar field ip besides the 
metric 17^^ (hence the index 1 on Ti{ao, /?o))- Second, it contains only positive- 
energy excitations (while one combination of the two scalar fields of T2(/3', /3") 
carried negative-energy waves). Third, it is the minimal way to parametrize 
the huge class of tensor-mono-scalar theories with a 'coupling function' a{(p) 
satisfying some very general requirements (see below). 

Let us now motivate the use of tensor-scalar theories of gravity as alternatives 
to general relativity. 

5.3 Tensor-scalar theories of gravity 

Let us start by recalling (essentially from 35j) why tensor- (mono)-scalar theories 
define a natural class of alternatives to GR. First, and foremost, the existence 
of scalar partners to the graviton is a simple theoretical possibility which has 
surfaced many times in the development of unified theories, from Kaluza-Klein 
to superstring theory. Second, they are general enough to describe many inter- 
esting deviations from GR (both in weak- field and in strong field conditions), 
but simple enough to allow one to work out their predictions in full detail. 
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Let us therefore consider a general tensor-scalar action involving a metric 
dtiiy (with signature 'mostly plus'), a scalar field $, and some matter variables 
iprn (including gauge bosons): 



167rG* 



.... 

For simplicity, we assume here that the weak equivalence principle is satisfied, 
i.e., that the matter variables V'm are all coupled to the same 'physical metricr i 
g^i,. The general model (j40p involves three arbitrary functions: a function F{^) 
coupling the scalar $ to the Ricci scalar of 5^1^, 7? = R{g^^), a function Z{^) 
renormalizing the kinetic term of <i>, and a potential function [/($). As we 
have the freedom of arbitrary redefinitions of the scalar field, $—!■$' = /(<&), 
only two functions among F, Z and U are independent. It is often convenient 
to rewrite ([3D]) in a ccinoTiicul form, obtained by redefining both $ and 
according to 

9;.^Fi<P)~g,,, (41) 



ip = ± 



3 f '^($) 1 Z($) i^/^ 

4 i^2($) + 2 F($)_ 



(42) 



This yields 



S = J ^ gV' [R* - 25r5M^9.</' - Viip)] + Sra [V™; A'(^) gU , 

(43) 

where i?* = R{g'^^), where the potential 

V{ip)^F-^{'S?)U{'S?), (44) 
and where the conformal coupling function A{lp) is given by 

A{^) ^ F-^'\<^) , (45) 



with '^{(f) obtained by inverting the integral 

The two arbitrary functions entering the canonical form (|43l) are: (i) the con- 
formal coupling function A{Lp), and (ii) the potential function V{if>). Note that 
the 'physical metric' g^i, (the one measured by laboratory clocks and rods) is 
conformally related to the 'Einstein metric' i?*^, being given by g^^ = A^{if) 5*^. 
The canonical representation is technically useful because it decouples the two 



Actually, most unified models suggest that there are violations of the weak equivalence 
principle. However, the study of general string-inspired tensor-scalar models [71] has found 
that the composition-dependent effects would be negligible in the gravitational physics of 
neutron stars that we consider here. The experimental limits on tests of the equivalence 
principle would, however, bring a strong additional constraint of order 10~^ Oq ~ Aa/a < 
10~^^. As this constraint is strongly model-dependent, we will not use it in our exclusion 
plots below. One should, however, keep in mind that a limit on the scalar coupling strength 
Op of order Qq < 10~^ | 71l I72| is likely to exist in many, physically-motivated, tensor-scalar 
models. 
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irreducible propagating excitations: the spin-0 excitations are described by tp, 
while the pure spin-2 excitations are described by the Einstein metric 17*^ (with 
kinetic term the usual Einstein-Hilbert action cx R{g'^^))- 

In many technical developments it is useful to work with the logarithmic 
coupling function a{ip) such that: 

a{(p) = \nA{ip) ; A{ip) = e"^'^) . (46) 

In the case of the general model (|40l) this logarithmic^ coupling function is 
given by 

a(^) = -ilnF(<i>), 

where 3>((^) must be obtained from (j42|) . 

In the following, we shall assume that the potential V{(p) is a slowly varying 
function of ip which, in the domain of variation we shall explore, is roughly 
equivalent to a very small mass term V{ip) 2 ^^^{ip — po)'^ with of cosmo- 
logical order of magnitude = 0{Hq), or, at least, with a range = 
much larger than the typical length scales that we shall consider (such as the 
size of the binary orbit, or the size of the Galaxy when considering violations of 
the strong equivalence principle) . Under this assumptiorP^ the potential func- 
tion V{ip) will only serve the role of fixing the value of ip far from the system 
(to p{r — 00) — po)i E^nd its effect on the propagation of p within the system 
will be negligible. In the end, the tensor-scalar phenomenology that we shall 
explore only depends on one function: the coupling function a{p). 

Let us consider some examples to see what kind of coupling functions might 
naturally arise. First, the simplest case is the Jordan-Ficrz-Brans-Dicke action, 
which is of the general type (HOI) with 

f($) = $ (47) 

Z($) = WBD$^\ (48) 

where wbd is an arbitrary constant. Using Eqs. (|42|) . (|45|) above, one finds that 
— 2aoip = hi^ and that the (logarithmic) coupling function is simply 

a{p) — uqp + const. , (49) 

where Uq = =F(2ci;BD + 3)"^^^, depending on the sign chosen in Eq. (|42| . Inde- 
pendently of this sign, one has the link 

(50) 



2 UJBD + 3 



^^As we shall mostly work with a{ip) below, we shall henceforth drop the adjective 'loga- 
rithmic'. 

^^Note, however, that, as was recently explored in |73l 1741 175) . a sufficiently fast varying 
potential V{ip) can change the tensor-scalar phenomenology by endowing ip with a mass term 
= i d^V/dip^ which strongly depends on the local value of 1/3 and, thereby can get large 
in sufficiently dense environments. 



20 



Note that 2uj-bu + 3 must be positive for the spin-0 excitations to have the 
correct (non ghost) sign. 

Let us now discuss the often considered case of a massive scalar field having 
a nonminimal coupling to curvature 



16nG 

This is of the form (gOl) with 

i^($) ^ 1 + ^$2 



(51) 



Z($) = 1, [/($) =m|$2. 



(52) 



The case ^ = — ^ is usually referred to as that of 'conformal coupling'. With the 
variables (HI]) the theory is ghost-free only if 2 (1 + ^$^^2 (^^^pj^^y ^ i ^ ^(i + 
6 ^) is everywhere positive. If we do not wish to restrict the initial values of 
we must have ^(1 + 6^) > 0. Introducing then the notation x = \/ C(l + 6 ^), 
we get the following link between $ and ip: 



2V2LP = I In 
+ 76 In 



1 + 2 X $ ( VT+X^ + ^ $ 



1 - 2\/6^$ 



(53) 



For small values of $, this yields Lp = 
coupling functions are given by 



m| $2 
1 + ^$2 



1 + ^$2 
0($^). The potential and the 

(54) 



ln(l + ^$2 



(55) 



These functions have singularities when 1 + ^$2 vanishes. If we do not wish 
to restrict the initial value of $ we must assume ^ > (which then implies our 
previous assumption C(l + 6(^) > 0). Then there is a one-to-one relation between 
$ and ip over the entire real line. Small values of <1? correspond to small values 
of and to a coupling function 



a{tf) 



(56) 



On the other hand, large values of |<I>| correspond to large values of |(^|, and to 
a coupling function of the asymptotic form 



a{ip) 



-72^1^1 
X 



const. 



(57) 



The potential V{ip) has a minimum at = 0, as well as other minima at 
ip ±cxD. If we assume, for instance, that m\ and the cosmological dynamics 
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are such that the cosmological value of ip is currently attracted towards zero, 
the value of at large distances from the local gravitating systems we shall 
consider will he ipo <^1. 

As a final example of a possible tensor-scalar gravity theory, let us discuss 
the string-motivated dilaton-runaway scenario considered in [76j . The starting 
action (a functional of 5^1^ and was taken of the general form 



and it was assumed that all the functions Bi{^) have a regular asymptotic be- 
havior when $ +00 of the form Bi{^) — Ci + 0{e^'^). Under this assumption 
the early cosmological evolution can push $ towards +00 (hence the name 'run- 
away dilaton'). In the canonical, 'Einstein frame' representation (j43p . one has, 
for large values of $ ~ c (/9, where c is a numerical constant, and the coupling 
function to hadronic matter is given by 



where 63 is the one-loop rational coefficient entering the renormalization-group 
running of the gauge field coupling g^. This finally yields a coupling function 
of the approximate form (for large values of tp): 



where the dimensionless constants k and c are both expected to be of order unity. 
[The constant c must be positive, but the sign of k is not a priori restricted.] 

Summarizing: the JFBD model yields a coupling function which is a linear 
function of f, Eq. (|49)l . a nonminimally coupled scalar yields a coupling function 
which interpolates between a quadratic function of (p, Eq. (j56p . and a linear one, 
Eq. ([57]) , and the dilaton-runaway scenario of Ref. [75] yields a coupling function 
of a decaying exponential type. 

5.4 The role of the couphng function a{Lp)\ definition of the 
two-dimensional space of tensor-scalar gravity theo- 



Let us now discuss how the coupling function a{p) enters the observable pre- 
dictions of tensor-scalar gravity at the first post-Newtonian (IPN) level, i.e., 
in the weak-field conditions appropriate to solar-system tests. It was shown in 
previous work that, if one uses appropriate units in the asymptotic region far 
from the system, namely units such that the asymptotic value a((/3o) of a{^p) 
vanished, all observable quantities at the IPN level depend only on the values 

^^In these units the Einstein metric g"^^ and the physical metric asymptotically coin- 




e"('^) cx Aqcd(^) (X B-^/\^) exp[-8^2 ^-1 ^^(^^j 



a{Lp) ~ke + const. , 



ries Ti(ao,/5o) 



cide. 
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of the first two derivatives of the a((^) at (p — (po. More precisely, if one defines 

, , da(w) ^, , da(Lp) d'^ a{w) 
a = ■ /? = — ^ = — ^ , 58 

dip dip ap)^ 

and denotes by ao = a{ipQ\ Pq = P{ipo) their asymptotic values, one finds 
(see, e.g., [33i) that the effective gravitational constant between two bodies (as 
measured by a Cavendish experiment) is given by 

G^G,{l + al), (59) 

while, among the PPN parameters, only the two basic Eddington ones, 7 = 
^PPN _ -|^^ g^j^j p ^ /jPPN — i (]^Q not vanish, and are given by 



7-7^™-l = -2Tf^, (60) 
5 _ flPPN ^ 1 "0 Q^o fai^ 



The structure of the results (pOj) and (pTjl can be transparently expressed by 
means of simple (Feynman-like) diagrams (see, e.g., [77] )■ Eqs. ([55)1 and ((5(1)) 
correspond to diagrams where the interaction between two worldlines (repre- 
senting two massive bodies) is mediated by the sum of the exchange of one 
graviton and one scalar particle. The scalar couples to matter with strength 
ao \/G*. The exchange of a scalar excitation then leads to a term oc a§. On 
the other hand, Eq. (I6ip corresponds to a nonlinear interaction between three 
worldlines involving: (i) the 'generation' of a scalar excitation on a first world- 
line (factor ao), (ii) a nonlinear vertex on a second worldline associated to the 
quadratic piece of a{<p) (aquad(v') — 5 Pai'-P — "y^o)^; so that one gets a factor /3o), 
and (iii) the final 'absorption' of a scalar excitation on a third worldline (second 
factor ao)- 

Eqs. ([60| and ([6T|l can be summarized by saying that the first two coefficients 
in the Taylor expansion of the coupling function &(</?) around ip = pQ (after 
setting a[po) — 0) 

a{ip) ^ aa{ip - Pa) +'^(}o{'^ - Pof ^ (62) 

suffice to determine the quasi-stationary, weak- field (IPN) predictions of any 
tensor-scalar theory. In other words, the solar-system tests only explore the 
'osculating approximation' (|62p (slope and local curvature) to the function a{p). 
Note that GR corresponds to a vanishing coupling function a{p) = (so that 
Q^o = /3o = • • • = 0), the JFBD model corresponds to keeping only the first term 
on the R.H.S. of ((62|) . while, for instance, the nonminimally coupled scalar field 
(with asymptotic value (po 1) does indeed lead to nonzero values for both ao 
and /3o, namely 

ao~-2e(^o; /3o--2e. (63) 
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Finally the dilaton-runaway scenario considered above leads also to non zero 
values for both ao and /3o, namely 

ao ~ - /c c e-"'^» ; (3o-+kc^ e'^'f' , (64) 

for a largish value of ipQ. Note that the dilaton-runaway model naturally predicts 
that ao ^ 1, and that /3o is of the same order of magnitude as ao : /3o — — cao 
with c being (positive and) of order unity. The interesting outcome is that such 
a model is well approximated by the usual JFBD model (with /3o = 0). This 
shows that a JFBD-like theory could come out from a model which is initially 
quite different from the usual exact JFBD theory. 

As we shall discuss in detail below, solar-system tests constrain ag and a§ |/3o| 
to be both small. This immediately implies that |ao| must be small, i.e., that 
the scalar field is linearly weakly coupled to matter. On the other hand, the 
quadratic coupling parameter /3o is not directly constrained. Both its magnitude 
and its sign can be more or less arbitrary. Note that there are no a priori sign 
restrictions on /3o. The conformal factor A^{lp) = exp(2a((/?)) entering Eq. I|43p 
had to be positive, but this leads to no restrictions on the sign of a{ip) and of 
its various derivatively. For instance, in the nonminimally coupled scalar field 
case, it seemed more natural to require f > 0, which leads to a negative /3o in 
view of Eq. 

Let us summarize the results above: (i) the most general tensor-scalar the- 
orjf^is described by one arbitrary function a((p); and (ii) weak- field tests depend 
only on the first two terms, parametrized by ao and /3o, in the Taylor expansion 
(j62p of a{(p) around its asymptotic value Lp^. 

From this follows a rather natural way to define a simple mini space of tensor- 
scalar theories. It suffices to consider the two-dimensional space of theories, say 
Ti(ao,/9o), defined by the coupling function which is a quadratic polynomial in 
[MIES], say 

aaoAX'P) = ao{ip - ipo) + -^M^ - ^af ■ (65) 

As indicated, this class of theories depends only on two parameters: ao and 
/3o- The asymptotic value (/Jq of does not count as a third parameter (when 
using the form (|65|) ) because one can always work with the shifted field = 
if — ipo, with asymptotic value (po = and coupling function aag^/s^^ {(p) = ao(p + 
I /3o t^^- Moreover, as already said, the asymptotic value a((/?o) of a{ip) has also 
no physical meaning, because one can always use units such that it vanishes (as 
done in ([^ ). 



^^As explained above, we assume here the presence of a potential term V{ifi) to fix the 
asymptotic value ipo of ip. If the potential V{ip) is absent (or negligible), the 'attractor 
mechanism' of Refs. 1781 1711 would attract to a minimum of the coupling function a{tp), 
thereby favoring a positive value of /Jq. 

-"^^Under the assumption that the potential Vi^ip) is a slowly-varying function of i/p, which 
modifies the propagation of i/s only on very large scales. 
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Note also that an alternative way to represent the same class of theories is 
to use a coupling function of the very simple form 

aM^) = ^/3(p^ (66) 

but to keep the asymptotic value as an independent parameter. This class 
of theories is clearly equivalent to Ti(aQ^ (3a), Eq. (j65p . with the dictionary: 

ao = /3<Po, Po = (3- 

5.5 Tensor-scalar gravity, strong-field effects, and binary- 
pulsar observables 

Having chosen some mini-space of gravity theories, we now wish to derive what 
predictions these theories make for the timing observables of binary pulsars. To 
do this we need to generalize the general relativistic treatment of the motion and 
timing of binary systems comprising strongly self-gravitating bodies summarized 
above. Let us recall that this treatment was based on a multi-chart method, 
using a matching between two separate problems: (i) the 'internal problem' 
considers each strongly self-gravitating body in a suitable approximately freely 
falling frame where the influence of its companion is small, and (ii) the 'external 
problem' where the two bodies are described as effective point masses which 
interact via the various fields they are coupled to. Let us first consider the 
internal problem, i.e., the description of a neutron star in an approximately 
freely falling frame where the infiuence of the companion is reduced to imposing 
some boundary conditions on the tensor and scalar fields with which it interacts 
[7l[8l[33l[34l[35]. The field equations of a general tensor-scalar theory, as derived 
from the canonical action (|43p (neglecting the effect of V{(p)) read 

R;, = 2 d,^ d,^ + Svr G, (t;, - i T*g;,^ , (67) 

□g. ^ = -4:TTG^a{ip)n, (68) 

where T^" = 2 c (g*)~^^^ 5Sra/5g*^^ denotes the material stress-energy tensor in 
'Einstein units', and a{Lp) the 95-derivative of the coupling function, see Eq. (|58[) . 
All tensorial operations in Eqs. (j67p and (|68[) are performed by using the Einstein 
metric g^^. 

Explicitly writing the field equations (p7|) and ([55)1 for a slowly rotating 
(stationary, axisymmetric) neutron star, labellecF"! A, leads to a coupled set of 
ordinary differential equations constraining the radial dependence of 17*^ and 
(f> [35t [79] . Imposing the boundary conditions g*^ rjf^,^, ip ^ ipa dX large 
radial distances, finally determines the crucial 'form factors' (in Einstein units) 
describing the effective coupling between the neutron star A and the fields to 

^'^We henceforth use the labels A and B for the (recycled) pulsar and its companion, instead 
of the labels a and b used above. We henceforth use the label a to denote the asymptotic 
value of some quantity (at large radial distances within the local frame, or Xg, of the 
considered neutron star A or B). 
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which it is sensitive: total mass ruAifa), total scalar charge ojAifa), and inertia 
moment J^(</3q). As indicated, these quantities are functions of the asymptotic 
value ipa of if felt by the considered neutron staiF^. They satisfy the relation 
LVAifa) — ^dmAifa)/dipa- From them, one defines other quantities that play 
an important role in binary pulsar physics, notably 

(jJa _ dlnniA 

mA~ difia 



aAi^Pa) 



(69) 



as well as 



kA{^a) = 



daA 
d(pa ' 

d\nlA 

dip a 



(70) 



(71) 



The quantity a a, Eq- ((69)) . plays a crucial role. It measures the effective coupling 
strength between the neutron star and the ambient scalar field. If we formally let 
the self-gravity of the neutron A tend toward zero (i.e., if we consider a weakly 
self-gravitating object), the function aA(}Pa) becomes replaced by a{ipa) where 
a((^) = da{(p)/d(p is the coupling strength appearing in the R.H.S. of Eq. ([55]) . 
Roughly speaking, we can think of UAifa) as a (suitable defined) average value 
of the local coupling strength a{(p{r)) over the radial profile of the neutron star 
A. 



neutron star 



a- 



scalar charge 
IkaI 




baryonic mass 
mA/mg 



critical maximum maximum 

mass mass in GR mass 



Figure 2: Dependence upon the baryonic mass ttia of the coupling parameter 
ttA in the theory Ti{ao,l3o) with ao = —0.014, /3o = —6. Figure taken from 

m- 

It was pointed out in Refs. [SU [35] that the strong self-gravity of a neutron 
star can cause the effective coupling strength aAi'fia) to become of order unity, 

■^^This ipa is a combination of the cosmological background value ipo and of the scalar 
influence of the companion of the considered neutron star. It varies with the orbital period and 
is determined as part of the 'external problem' discussed below. Note that, strictly speaking, 
the label a (for asymptotic) should be indexed by the label of the considered neutron star: i.e. 
one should use a label (and a locally asymptotic value ) when considering the neutron 
star A, and a label ag (with a corresponding ^Pub) when considering the neutron star B. 
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even when its weak-field counterpart ao = Oi{'Pa) is extremely small (as is im- 
plied by solar-system tests that put strong constraints on the PPN combination 
7 = — 2 ttp/ (1-l-aQ)). This is illustrated, in the minimal context of the Ti[aQ, f3o) 
class of theories, in Figured) 

Note that when the baryonic mass of the neutron star is smaller than 
the critical mass fficr — 1.24 Mq the effective scalar coupling strength a a of 
the star is quite small (because it is proportional to its weak-field limit ao = 
a{(pa))- By contrast, when rfiA > rhcr, \o:a\ becomes of order unity, nearly 
independently of the externally imposed ao = aa = a{ipa)- This interesting 
non-perturhative behaviour was related in |34[|35| to a mechanism of spontaneous 
scalarization, akin to the well-known mechanism of spontaneous magnetization 
of ferromagnets. See also [51^ for a simple analytical description of the behaviour 
of aA- 

Let us also mention in passing that, in the case where ^4 is a black hole, the 
effective coupling strength a a actually vanishes [33]. This result is related to 
the impossibility of having (regular) 'scalar hair' on a black hole. 

We have sketched above the first part of the matching approach to the mo- 
tion and timing of strongly self-gravitating bodies: the 'internal problem'. It 
remains to describe the remaining 'external problem'. As already mentionned 
(and emphasized, in the present context, by Eardley [3 HI]), the most efficient 
way to describe the external problem is, instead of matching in detail the exter- 
nal fields (g*^ , if) to the fields generated by each body in its comoving frame, to 
'skeletonize' the bodies by point masses. Technically this means working with 
the action 

S = I ^ g'AR. - 2 gr dM 

- f mA{^{zA)){-9l,{zA)dz>Xdz'XY'\ (72) 

A •' 

where the function mAi'fi) in the last term on the R.H.S. is the function myi((^a) 
obtained above by solving the internal problem. Eq. ([72l) indicates that the ar- 
gument of this function is taken to be ipa = ip{zA), i.e., the value that the scalar 
field (as viewed in the external problem) takes at the location of the center of 
mass of body A. However, as body A is described, in the external problem, as a 
point mass this causes a technical difficulty: the externally determined field ip{x) 
becomes formally singular at the location of the point sources, so that ip{za) is a 
priori undefined. One can either deal with this problem by coming back to the 
physically well-defined matching approach (which shows that f{zA) should be 
replaced by (pa, the value of ip in an intermediate domain Ra <C r ^ \za — Zb\), 
or use the efficient technique of dimensional regularization. This means that the 
spacetime dimension D in Eq. ([7^ is first taken to have a complex value such 
that f{zA) is finite, before being analytically continued to its physical value 
D = 4. 

One then derives from the action ((72]) two important consequences for the 
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motion and timing of binary pulsars. First, one derives the Lagrangian de- 
scribing the relativistic interaction between N strongly self-gravitating bodies 
(including orbital ~ (w/c)^ effects, and neglecting ©(ti^/c^) ones) [TTl [71 [551 [55] . 
It is the sum of one-body, two-body and three-body terms. 

The one-body action has the usual form of the sum (over the label A) of the 
kinetic term of each point mass: 



9 1 ?! {vl? fl\ 
= -mAC^ + -mAv\ + -mA^^ + 0{^\. (73) 

Here, we use Einstein units, and the inertial mass entering Eq. (|73p is tua = 
mA{vo), where lpq is the asymptotic value of ip far away from the considered 
iV-body system. 

The two-body action is a sum over the pairs A,Boi& term L^^'^^ which 
differs from the GR-predicted 2-body Lagrangian in two ways: (i) the usual 
gravitational constant G appearing as an overall factor in L^ab'^^ must be re- 
placed by an effective (body-dependent) gravitational constant (in the appro- 
priate units mentioned above) given by 

GAi3 = G^l + a^as), (74) 

and (ii) the relativistic {Oiv^ / c^)) terms in L^^'^^ contain, in addition to those 
predicted by GR, new velocity-dependent terms of the form 

r-yr 2-body ,- sGABrUAmB {va-Vb)'^ 

S^L^B -ilAB) — , (75) 

with 

lAB = lAB - 1 = - 2 — . (76) 

1 + UAOtB 

In these expressions a a = ola{pq) = dhmiAiva) / dkpi) (see Eq. (|69p with ipa 
Vo)- 

Finally, the 3-body action is a sum over the pairs B, C and over A (with 
A ^ B, A ^ C, but the possibility of having B = C) of 

r 3-body _ ^ Gab Gac ^a ttib mg , . 

^ABC --(1 + 2/W c^rABTAC ^ ^ 



where 



nA ^ r,A 1 _ 1 (^A ac 

Pbc = Pbc ^ ^ - o 77~, Tf;—, T ' \'°) 

2 {l + aA aB){l + aA ac) 

with Pa = daA{'Po)/dipQ (see Eq. ([70)) with Lpa Vo)- 

When comparing the strong-field results ((74)) . ((76)) . (|78|) to their weak-field 
counterparts ([59]), (|60]) . (|6T]l one sees that the body-dependent quantity ua 
replaces the weak-field coupling strength uq in all quantities which are linked 
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to a scalar effect generated by body A. Note also that, in keeping with the 
'3-body' nature of Eq. ([77)) . the quantity (3^q — 1 is linked to scalar interactions 
which are generated in bodies B and C and which nonlinearly interact on body 
A. The notation used above has been chosen to emphasize that 'jab and 
are strong- field analogs of the usual Eddington parameters 7^^"^, so that 

7a_b and (3gQ are strong-field analogs of the 'post-Einstein' IPN parameters 7 
and P (which vanish in GR). Indeed the usual PPN results for the post-Einstein 
terms in the 0(l/c^) 2-body and 3-body Lagrangians are obtained by replacing 
in Eqs. ^ and ^ ^ab 7, Pbc P and Gab G. 

The non-perturbative strong-field effects discussed above show that the strong 
self-gravity of neutron stars can cause jab and to be significantly different 



from their GR values 7' 



GR _ 



1, 



GR _ 



1, in some scalar-tensor theories having 



— 1 cx ttQ and 



a small value of the basic coupling parameter ao (so that 7^^^ 
^PPN — I (X Pq are both small). For instance, Fig. O shows that it is possible 
to have aA ~ cub ± 0.6 which implies jab — 1 ~ — 0.53, i.e., a 50% deviation 
from GR! Even larger effects can arise in — 1 because of the large values 
that Pa = duA/dipQ can reach near the spontaneous scalarization transition 



Those possible strong-field modifications of the effective Eddington param- 
eters Jab, PbC' which parametrize the 'first post-Keplerian' (IPK) effects 
(i.e., the orbital effects ^ jt? smaller than those entailed by the Lagrangian 
vriA v\ + I X] Gab rriA mB/fAB), can then significantly modify the usual 

A A^B 

GR predictions relating the directly observable parametrized post-Keplerian 
(PPK) parameters to the values of the masses of the pulsar and its compan- 
ion. As worked out in Refs. [TTJ [3TJ [33l [35] one finds the following modified 
predictions for the PPK parameters fc = r and s: 



k*-^{mA,mB) = 



GAB{mA + rns) n 



2/3 



l-\aAaB _ Xa Pb a\ + Pa a% 
1 + aACeB 6 (1 + Q!A Q:_b)2 



(79) 



„th 



{mA,mB) = Gobtub , 



„th 



{mA,mB) 



n XA 



GAB{mA + tob) 71 



-1/3 



(80) 
(81) 



Here, the label A refers to the object which is timed ('the pulsar ri). the label B 
refers to its companion, xa = a^sini/c denotes the projected semi- major axis 
of the orbit of A (in light seconds), Xa = "t-a/ {mA-'r'niB) and Xb = ms / {'niA + 
ruB) = 1 — Xa the mass ratios, n = 2t:/Pi, the orbital frequency and Gqb = 
G'*(l -|- Qfocts) the effective gravitational constant measuring the interaction 



■^^In the double binary pulsar, both the first discovered pulsar and its companion are pulsars. 
However, the companion _B is a non recycled, slow pulsar whose motion is well described by 
Keplerian parameters only. 
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between B and a test object (namely electromagnetic waves on their way from 
the pulsar toward the Earth). In addition one must replace the unknown bare 
Newtonian G« by its expression in terms of the one measured in Cavendish 
experiments, i.e., G* = G/(l + a^) as deduced from Eq. ([55)1 . 

The modified theoretical prediction for the PPK parameter 7 entering the 
'Einstein time delay' A^, Eq. is more complicated to derive because one 

must take into account the modulation of the proper spin period of the pulsar 
caused by the variation of its moment of inertia I a under the (scalar) influence 
of its companion [TTl [7l [35] . This leads to 

th/ N e Xb f GABimA + mB)n^^^^ 

7 (mA,mB) - ' 



n 1 + a A (Xb \ 
[XBil + aAae) + l + kAaB], (82) 

where fc^(</9o) = ~dhilA{^o)/dvo (see Eq. ([7T|) with ipa ifio). Numerical 
studies |3S] show that kA can take quite large values. Actually, the quantity 
kACtB entering ()82[) blows up near the scalarization transition when ao — > 
(keeping Po < fixed). In other words a theory which is closer to GR in weak- 
field conditions predicts larger deviations in the strong-field regime. 

The structure dependence of the effective gravitational constant Gab , Eq. ([74] 
has also the consequence that the object A does not fall in the same way 
as B in the gravitational field of the Galaxy. As most of the mass of the 
Galaxy is made of non strongly-self-gravitating bodies, A will fall toward the 
Galaxy with an acceleration oc Gaq, while B will fall with an acceleration 
oc Gbq- Here, as above, Gao — Gqa — G*(I 4- aoUA) is the effective gravi- 
tational constant between A and any weakly self-gravitating body. As pointed 
out in Ref. [43] this possible violation of the universality of free fall of self- 
gravitating bodies can be constrained by using observational data on the class 
of small-eccentricity long-orbital-period binary pulsars. More precisely, the 
quantity which can be observationally constrained is not exactly the violation 
A AS = {Gqa — Gqb) I G — ( I -f ag ) ~ ^ (ao ttA — cto ckb ) of the strong equivalence 
principle [which simplifies to A^o = {Gqa — G)/G = (1 + a'^)^^{aQaA — oiq) 
in the case of observational relevance where one neglects the self-gravity of the 
white-dwarf companion] but rathe@ [55] 

^effective ^ ( 1 AB - {X A PaA + X B Pb b) + 



3 

(I -I- aA as)"^^^(f + aly^{aQ ua - ao as) ■ (83) 

Here, the index B {— white-dwarf companion) can be replaced by (weakly self- 
gravitating body) so that, for instance, jab — 7ao = 1 — 2 a a Q!o/(1-I- a a olq) = 
(1 — aA Q!o)/(l + aA ao), as deduced from Eq. ([76]) . 



■^•^This refinement is given here for pedagogical completeness. However, in practice, the 
lowest-order result A ~ (1 + a^^~^(ci^^ — a^) ~ oljs^ — is accurate enough. 
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It remains to discuss the possible strong-field modifications of the theoretical 
prediction for the orbital period derivative Ph = P^^{mA,'mB)- This is obtained 
by deriving from the effective action ([7^ the energy lost by the binary system 
in the form of fluxes of spin-2 and spin-0 waves at infinity. The needed results 
in a generic tensor-scalar theory were derived in Refs. [531 132] (in addition one 
must take into account the tensor-scalar modification of the additional 'varying- 
Doppler' contribution to the observed Pb due to the Galactic acceleration [38]). 
The final result for Pf, is of the form 

pth ('^ ^ ^ _ pinonopolc . Adipolc , Aquadrupolo . Aquadrupolc 

-Tf, (mA,mB) — r-i,^ + ^btp + ^b^p ^bg' 

+ ps^i^s"= + ,5*'^Ff'"^"% (84) 

where, for instance, p^^°P°^° jg (heuristicalljQ) related to the monopolar flux 
of spin-0 waves at infinity. The term p^^^'^^^'£'°^^ corresponds to the usual 
quadrupolar flux of spin-2 waves at inflnity. It reads: 



Aquadrupolc/ x 1927r TO^ TOb f^^. 

Pbg' {mA,mB) = -777- TT 7 yl (85) 

(GAB{mA + mB)nV'^ 1 + 73 eV24 + 37 eV96 



V c3 y (l-e2)7/2 

with Gab — G,(l -t- as) — G{1 + aA a_B)/(l + Oq), where G, is the 'bare' 
gravitational constant appearing in the action, while G is the gravitational con- 
stant measured in Cavendish experiments. The flux (j85[) is the only one which 
survives in GR (although without any a^^-related modifications). Among the 
several other contributions which arise in tensor-scalar theories, let us only write 
down the explicit expression of the contribution to (j84p coming from the dipolar 
flux of scalar waves. Indeed, this contribution is, in most cases, the dominant 
one [7] because it scales as {v/cf , while the monopolar and quadrupolar con- 
tributions scale as {v/c)^. It reads 

pdipoio/ X „ G^mAmBTi 1 + 6^/2 2 /qrn 

PbJ {mA,mB) = ~2Tr— ; ^-ttj (a^ - as) • (86) 



c^{mA + tub) (1 — e^)^/^ 



bip 

Note that the dipolar effect ([SB)) vanishes when a a = ctB- Indeed, a binary 
system made of two identical objects {A = B) cannot select a preferred direction 
for a dipole vector, and cannot therefore emit any dipolar radiation. This also 
implies that double neutron star systems (which tend to have ttla ~ rriB ^ 
1.35 Mq) will be rather poor emitters of dipolar radiation (though ([55]) still tends 
to dominate over the other terms in ([M)) . because of the remaining difference 
(jua — rnB) / {ttla + rriB) 7^ 0). By contrast, very dissymmetric systems such 



Contrary to the GR case where a lot of effort was spent to show how the observed Pj, 
was directly related to the GR predictions for the (t)/c)^-accurate orbital equations of motion 
of a binary system 9 , we use here the indirect and less rigorous argument that the energy 
flux at infinity should be balanced by a corresponding decrease of the mechanical energy of 
the binary system. 
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as a neutron-star and a white-dwarf (or a neutron-star and a black hole) will 
be very efficient emitters of dipolar radiation, and will potentially lead to very 
strong constraints on tensor-scalar theories. See below. 

5.6 Theory-space analyses of binary pulsar data 

Having reviewed the theoretical results needed to discuss the predictions of 
alternative gravity theories, let us end by summarizing the results of various 
theory-space analyses of binary pulsar data. 

Let us first recall what are the best, current solar-system limits on the two 
IPN 'post-Einstein' parameters 7 = 7^^^ — 1 and (3 = (3^^^ — 1. They are: 

7= (2.1±2.3) X 10-^ (87) 

from frequency shift measurements made with the Cassini spacecraft [81) , which 
supersedes the constraint 

7 = (-1.7 ±4.5) X 10"'' (88) 
from VLBI measurements [82], 

127-^1 <3xl0"^ (89) 
from Mercury's perihelion shift [66|, I83j , and 

4/3-7= (4.4± 4.5) X 10"'', (90) 

from Lunar laser ranging measurements [84j . 

Concerning binary pulsar data, we can make use of the published measure- 
ments of various Keplerian and post-Keplerian timing parameters in the binary 
pulsars: PSR 1913-^16 [37], PSR B1534-hl2 [41], PSR J1141-6545 47J and 
PSR J0737-3039A-f-B f3l, gi [49]. In addition, we can useEl the recently up- 
dated limit on the parameter A measuring a possible violation of the strong 
equivalence principle (SEP), namely |A| < 5.5 x 10""^ at the 95% confidence 
level [H]. 

This ensemble of solar-system and binary-pulsar data can then be analyzed 
within any given parametrized theoretical framework. For instance, one might 
work within 

(i) the 4-parameter framework To (7, P; e, Q [70] which defines the 2PN exten- 
sion of the original (Eddington) PPN framework To(7,/3); or 

(ii) the 2-parameter class of tensor- mo no-scalar theories ri(ao,/3o) [34]; or 

■^^There is, however, a caveat in tlie theoretical use one can make of the phenomenological 
Umits on A. Indeed, in the small-eccentricity long-orbital-period binary pulsar systems used 
to constrain A one does not have access to enough PK parameters to measure the pulsar mass 
ruA directly. As the theoretical expression of A ~ oq CfA — ctp depends on rriA (through a a), 
one needs to assume some fiducial value of niA (say tua ~ 1.35 Mq). 
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(iii) the 2-parameter class of tensor-bi-scalar theories T2{/3' , f3") [53]. 

Here, the index on Tq{j, (3; e, C) is a reminder of the fact that this framework 
is not a family of specific theories (it contains zero explicit dynamical fields), 
but is a parametrization of 2PN deviations from GR. As a consequence, its use 
for analyzing binary pulsar data is somewhat ill-defined because one needs to 
truncate the various timing observables (which are functions of the compactness 
of the two bodies A and B, say P^^ — /{ca^cb)) at the 2PN order (i.e. es- 
sentially at the quadratic order in ca and/or cb)- For some observables (or for 
product of observables) there might be several ways of defining this truncation. 
In spite of this slight inconvenience, the use of the To(7,/3;e, C) framework is 
conceptually useful because it shows very clearly why and how binary-pulsar 
data can probe the behaviour of gravitational theories beyond the usual IPN 
regime probed by solar-system tests. 

For instance, the parameter A^i = m^'^^/m^'^''* — 1 measuring the strong 
equivalence principle (SEP) violation in a neutron star has, within the ro(7, /?; e, C) 
framework, a 2PN-order expansion of the form [33l [70] 

Aa = -^(4^-7)ca+(|+C + 0(^)) (91) 

where ca = -2 ~ ^ {U)a, bA ^ ^ {U^)a ^ B c\, with B ~ 1.026 and 

CA — kmA/MQ with k ~ 0.21. The general result (|9ip is compatible with the 
result quoted in subsection 15.21 within the context of the theory T2{/3' , j3") when 
taking into account the fact that, within T2(/3',/3"), one has /3 = 7 = 0, e = /3' 
and C = [and that /?" parametrizes some effects beyond the 2PN level]. 

On the example of Eq. ([9T|l one sees that, after having used solar-system tests 
to constrain the first contribution on the RHS to a very small value, one can 
use binary-pulsar tests of the SEP to set a significant limit on the combination 
i e -I- C of 2PN parameters. Other pulsar data then yield significant limits on 
other combinations of the two 2PN parameters e and C- The final conclusion 
is that binary-pulsar data allow one to set significant limits (around or better 
than the 1% level) on the possible 2PN deviations from GR (in contrast to 
solar-system tests which are unable to yield any limit on e and Q [70]. For a 
recent update of the limits on e and which makes use of recent pulsar data 
see [51] . 

Let us now briefiy discuss the use of mini-space of theories, such as Ti(ao, /3o) 
or T2(/3',/3"), for analyzing solar-system and binary-pulsar data. The basic 
methodology is to compute, for each given theory (e.g. for each given values of 
ao and (3q if one chooses to work in the Ti(ao,/?o) theory space) a goodness- 
of-fit statistics x^(cio,/3o) measuring the quality of the agreement between the 
experimental data and the considered theory. For instance, when considering 
the timing data of a particular pulsar, for which one has measured several PK 
parameters pi {i = 1, . . . ,n) with some standard deviations cr°^*', one defines, 
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for this pulsar 

n 

x'(ao,/3o)= min Y.{of;)-\pf''°'^a^, P^-mA,mB) - ptl\ (92) 

1=1 

where 'min' denotes the resuh of minimizing over the unknown masses niA, ms 
and where p*'^''™-^(ao, /?o; "^s) denotes the theoretical prediction (within 
Ti(ao,/3o)) for the PK observable pi (given also the observed values of the 
Keplerian parameters). 

The goodness-of-fit quantity x^(ao, /3o) will reach its minimum x^jj^ for some 
values, say a™'", of ao and /3o- Then, one focusses, for each pulsar, on the 
level contours of the function 

Ax'(ao,/5o) = x'(ao,/3o) -Xmin- (93) 

Each choice of level contour (e.g. Ax^ 1 or Ax^ = 2.3) defines a certain 
region in theory space, which contains, with a certain corresponding 'confidence 
level', the 'correct' theory of gravity (if it belongs to the considered mini-space of 
theories). When combining together several independent data sets (e.g. solar- 
system data, and different pulsar data) we can define a total goodness-of-fit 
statistics Xtot (o^o i /^o ) i by adding together the various individual X^ (o^o j /^o ) ■ This 
leads to a corresponding combined contour A Xtot(Q^Oi /3o)- 

Let us end by briefly summarizing the results of the theory-space approach 
to relativistic gravity tests. For detailed discussions the reader should consult 
Refs. [331 SOI 133 [50] , and especially the recent update [5T] which uses the 
latest binary-pulsar data. 

Regarding the two-parameter class of tensor-bi-scalar theories T2(/3', (3") the 
recent analysis f5T| has shown that the A x^iP' t P") corresponding to the double 
binary pulsar PSR J0737— 3039 was defining quite a small elliptical allowed 
region in the (/?', /3") plane. By contrast the other pulsar data define much wider 
allowed regions, while the strong equivalence principle tests define (in view of 
the theoretical result A ~ 1 -f- i B(3'{c\ — c^)) a thin, but infinitely long, strip 
l/3'l < est. in the (/?', (3") plane. This highlights the power of the double binary 
pulsar in probing certain specific strong-field deviations from GR. 

Contrary to the T2(/3',/?") tensor-bi-scalar theories, which were constructed 
to have exactly the same first post-Newtonian limit as GF@ (so that solar- 
system tests put no constraints on (3' and /3"), the class of tensor- mono-scalar 
theories Ti(aoi/9o) is such that its parameters ao and /?o parametrize both the 
weak-field IPN regime (see Eqs. ([60]) and ([6T|) above) and the strong-field regime 
(which plays an important role in compact binaries) . This means that each class 
of solar-system data (see Eqs. ([87| -([90 |) above) will define, via a corresponding 
goodness-of-fit statistics of the type, say 

^2 ri \ / Cassini\— 2 /;r,thcory / n \ -Cassini\2 

XCassini("0,P0) = (^T- ) (7 '{olO,Po)-l ) 

■^^However, this could be achieved only at the cost of allowing some combination of the two 
scalar fields to carry a negative energy flux. 
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a certain allowed regioiFI in the (ao,/9o) plane. As a consequence, the analysis 
in the framework of the Ti{ao,/3o) space of theories allows one to compare 
and contrast the probing powers of solar-system tests versus binary-pulsar tests 
(while comparing also solar-system tests among themselves and binary-pulsar 
ones among themselves). The result of the recent analysis [51] is shown in 
Figure [31 




Figure 3: Solar-system and binary-pulsar constraints on the two-parameter fam- 
ily of tcnsor-mono-scalar theories Ti(aQ, (3q). Figure taken from [51j . 

In Fig. [3l the various solar-system constraints ([87)) - ((90)) are concentrated 
around the horizontal Po axis. In particular, the high-precision Cassini con- 
straint is the lower small grey strip. The various pulsar constraints arc labelled 
by the name of the pulsar, except for the strong equivalence principle constraint 
which is labelled SEP. Note that General Relativity corresponds to the origin 
of the (ao,/3o) plane, and is compatible with all existing tests. 

The global constraint obtained by combining all the pulsar tests would, to a 
good accuracy, be obtained by intersecting the various pulsar-allowed regions. 
One can then see on Fig. [3] that it would be comparable to the pre-Cassini 
solar-system constraints and that its boundaries would be defined successively 
(starting from the left) by 1913-1-16, 1141-6545, 0737-3039, 1913-^16 again 
and 1141—6545 again. 

A first conclusion is therefore that, at the quantitative level, binary-pulsar 
tests constrain tensor-scalar gravity theories as strongly as most solar-system 

■^^Actually, in the case of the Cassini data, as it is quite plausible that the positive value of 
the published central value ^Casssim _ ^2.1 X 10~^ is due to unsubtracted systematic effects, 
we use (jC^assini _ 2 3 X 10~^ but ^Cassim _ Q Otherwise, we would get unreasonably strong 
1(T limits on Oq because tensor-scalar theories predict that 7 must be negative, see Eqs. II6OII 
and |[6T||. 
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tests (excluding the exceptionally accurate Cassini result which constrains Oq 
to be smaller than 1.15 x 10~^, i.e. |ao| < 3.4 x 10"'^). A second conclusion 
is obtained by comparing the behaviour of the solar-system exclusion plots and 
of the binary-pulsar ones around the negative f3o axis. One sees that binary- 
pulsar tests exclude a whole domain of the theory space (located on the left 
of (3o < —4) which is compatible with all solar-system experiments (even when 
including the very tight Cassini constraint). This remarkable qualitative feature 
of pulsar tests is a direct consequence of the existence of (non-perturbative) 
strong-field effects which start developing when the product — /3o ca (with ca 
denoting, as above, the compactness of the pulsar) becomes of order unity. 

6 Conclusion 

In conclusion, we hope to have convinced the reader of the superb opportunities 
that binary pulsar data offer for testing gravity theories. In particular, they 
have been able to go qualitatively beyond solar-system experiments in probing 
two physically important regimes of relativistic gravity: the radiative regime 
and the strong-field one. Up to now. General Relativity has passed with flying 
colours all the radiative and strong-field tests provided by pulsar data. However, 
it is important to continue testing General Relativity in all its aspects (weak- 
field, radiative and strong-field). Indeed, history has taught us that physical 
theories have a limited range of validity, and that it is quite difficult to predict 
in which regime a theory will cease to be an accurate description of nature. Let 
us look forward to new results, and possibly interesting surprises, from binary 
pulsar data. 
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